On the Spin-Hall effect in a dirty Rashba semiconductor 
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The impact of the boundary condition on the solution of the spin diffusion equations in a dirty 
Rashba semiconductor is reinvestigated. To this end first diffusion equations, which take into account 
the coupling between spin and charge, are derived. The investigation of the solution to the spin- 
diffusion equation near the boundary shows that the spin-current tensor vanishes discontinuously 
at the boundary. As a consequence, there is no spin-Hall effect in a dirty Rashba semiconductor. 



Spintronics is a rapidly developing field, which has the 
objective to learn more about the interaction of the elec- 
tron spin with its solid state environment. It is stim- 
ulated by the fact that the electron spin can also be 
manipulated electrically and thus be utilized in electri- 
cal devices. The interaction, which provides access to 
electric spin manipulation is the spin-orbit interaction, 
which is realized by the Rashba interaction in the sim- 
plest situationyj. 

The injection of spins into nonmagnetic semiconduc- 
tors, however, has proven to be a difficult task. There- 
fore, much attention is paid to the investigation of effects, 
which permit the electrical creation of non-vanishing 
magnetizations in the sample. Already the first inves- 
tigations on the impact of a lateral electric field on a 
two-dimensional Rashba semiconductor have shown that 
such effects in principle exist 0- According to Ref. 0, 
e.g., the lateral field leads to a homogeneous magnetiza- 
tion in the two-dimensional plane. Similar results have 
been obtained in the Refs.0 and recently. 

An effect, which might be useful in creating non- 
vanishing magnetizations, is the spin-Hall effect. The 
latter manifests itself in a non-vanishing magnetization 
at the sample boundaries in the presence of an electric 
field. In that it differs from the spin accumulation, which 
is present in the bulk. Whereas there seems to be some 
agreement over the magnitude of the spin accumulation 
in a Rashba semiconductor, the spin-Hall effect is still a 
rather controversial issue. On the one hand an intrinsic 
spin-Hall effect has been suggested in the Refs.|3] and jf| 
and investigated further in a number of papers. These au- 
thors argue, that the expection value of the spin-current 
operator, defined by the equation 



Jik = ^{Vi,crk}, 



(1) 



is non- vanishing in a clean Rashba semiconductor in the 
presence of an electric field F, and that therefore a mag- 
netization is created at the boundaries, (pi is the velocity 
operator, 07. are the Pauli matrices, i = 1,2, k = 1,2,3. 
The indices 1 and 2 and the labels x and y are used si- 
multaneously to characterize the vectors in the 2-d plane, 
the index 3 and the label z to characterize the perpen- 
dicular direction). The spin-Hall conductivity, defined 
by the ratio an = Jyz/F, is universal in clean systems, 



that is independent of the magnitude of the Rashba in- 
teraction. In dirty systems, in which the Rashba level 
splitting is small compared to the disorder energy, the 
spin-Hall conductivity vanishes 0, 0. Thus, it is ex- 
pected that there is no spin-Hall effect in such systems. 
On the other hand spin-Hall effect has been discussed in 
Ref. 0| m the hopping regime. Doing so, it has been 
assumed that the Rashba level splitting can be ignored. 
Thus, the system of Ref. 0] is always in the dirty limit. 

When investigating the differences between the ap- 
proaches of the Refs. 0-0 and the approach of Ref. 
[lfj we note the following: whereas a diffusion equation 
for the spin density has been derived in Ref.0] the con- 
clusions of the Refs. 0-0 are based only on the investiga- 
tion of the spin-current tensor l(T[l. However, it is unclear 
how the tensor Q is related to observable quantities^]. 
Therefore, there are no predictions of the magnitude of 
the magnetization at the sample boundaries. Moreover, 
the expection value of the spin-current is non- vanishing 
even in equilibrium ^]J- This fact raises the question 
whether the spin-current tensor is well defined at all . 
In Ref. [H}, by contrast, the spin-current tensor has been 
read off from the diffusion equation. To this end, the 
authors of Ref. lQj have taken the point of view, that all 
terms, which can be cast into the form of a derivative, 
contribute to the divergence of the spin-current tensor. 
Using this approach they calculate the solution to the 
spin-diffusion equation and make detailed predictions to 
the magnitude of the magnetization at the sample bound- 
aries. However, also this approach is not unproblematic. 
Since the diffusion equations yields only the divergence 
of the spin-current tensor, the tensor read off from the 
diffusion equation is only determined up to a curl and 
boundary currents had to be assumed to be non-existent 
at all. Moreover, the spin accumulation has entirely been 
ignored. 

Since there seems to be no rigorous method, which 
could be used to define the spin-current tensor unam- 
biguously, there seems to be no other way than to in- 
vestigate the consequences of the definition (JJJ and to 
compare them with experiments. Acording to Ref.|lfj|| 
the spin-current tensor determines the boundary con- 
ditions in the diffusive regime. Thus, it is natural to 
focus on boundary effects. To get further insight into 
the implications of the definition Q and into the ques- 
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tion, in which way the spin current tensor determines the 
boundary conditions, we investigate the magnetization in 
a dirty Rashba semiconductor in a half plane. To this end 
we derive spin-diffusion equations for a two dimensional 
system and solve these equations subjected to the bound- 
ary condition that the spin-current tensor vanishes at the 
boundary. 

Diffusion equations, which also take into account the 
coupling between spins and charges, have been first de- 
rived in the Refs.lj] and 0. Here, we use the same 
method as in Ref.[9j. The Hamiltonian is given by the 
equation 



H 



2m 



(er, TV x p) + V(x) 



(2) 



Here p is the momentum operator, TV = Ne z , and m 
is the effective mass. V(x) is a random potential with 
Gaussian distribution function, zero mean and standard 
deviation 



{{V(x)V(x'))) = 



2~KVT 



S(x - x 1 ). 



(3) 



v is the density of states per spin and r the single par- 
ticle relaxation time. The propagation of single particle 
excitations is described by the retarded and advanced 
Green's functions. These functions are the solution to 
the equation 

{{±i— + E)l - H)G R/A (x, x'\E, s) = lS{x - x'). (4) 

Here 1 is the 2x2 unit matrix, s is a fixed frequency 
and E is the total particle energy. For the calculation 
of their configuration averages we use the self-consistent 
Born approximation. In this approximation the Fourier- 
transform g R (k\E, s) of the configuration averaged re- 
tarded Green's functions G R (x — x'\E, s) is given by the 
equation 

g R (k\E,s) = g + (k\E,a)g-(k\E,a) 

xig-^kfos^-faNxp)), (5) 

where gg\k\E, s) = ihs/2 + E - E±(k) + ih/2T, E±(k) = 
h 2 k 2 /2m± \N\hk and g(k\E,s) = g+(k\E, s)\ N=0 . The 
Fourier transform g A (k\E, s) of the configuration aver- 
aged advanced Green's function G (x — x'\E,s) is ob- 
tained from Eq.(PJ by hermitian conjugation. 

The evolution of the particle density and the spin den- 
sity is described by a generalized diffusion equation. The 
densities are defined by the relationship 

f aa ,(x\E,t)=tr{tl>+(x,t)il) a ,{x,t)po(E)), (6) 

where tp£(x,t) and ip a (x,t) are creation and annihila- 
tion operators for particles with spin a, and po(E) is the 
initial density matrix. The particle density is given by 
the equation n(x\E,t) = trf(x\E,t), the spin density by 



the relationship Si(x\E,t) = tr(<Jif(x\E,t)). Note, that 
this definition differs from the conventional definition by 
a factor 2/h. The equation of motion for the quantity / 
takes the form 

J dx 1 T™i 2 (x,x 1 \E,s)f aia2 (x 1 \E,s) = f° a ,(x\E) (7) 

after Laplace transformation with respect to times (t — * 
s). A summation with respect to double indices has to 
be performed. The quantity f® a ,(x\E) is the initial con- 
dition. The kernel T in Eq.Q is given by the equation 

G R ai (x-x'\E,s)Gi 2a ,(x'-x\E,s)(8) 



2-evt 



in the ladder approximation. To simplify Eq.Q we 
use the hydrodynamics expansion. Doing so, we obtain 
the following diffusion equations after an inverse Laplace 
transformation 



C1TI 

— — DAn + Qt(N x V,S) 







7 £f 

— + SI-S - DAS -u.(N xV)x S 
at 

+ flrN xVn = 0. 



(9) 



(10) 



These equations agree with those of Ref. |9| . Here D = 
Er/m is the diffusion coefficient, lo s = AmD/Ti and 
n ik = nS ik (l + 5 i3 ), where = 4m 2 N 2 D/h 2 . 

To calculate the spin-current density we use the ex- 
pression 



J lk {x\E,s) 



2-KVT 



tr / dyG A (y -x\E,Q)Ji 



ik 



xG R (x-y\E,0)f(y\E, S ), (11) 



which is valid in the ladder approximation, provided we 
are interested in time scales which are large compared 
to T0|. To simplify the integral we use the same hydro- 
dynamic expansion as before. Doing so, we obtain the 
following expressions for components of the spin-current 
tensor 



Nlo s 

J xx — DX7 x Sx ~z S z 



J xy = —DV x S y + flrNn, 



(12) 



(13) 



Jxz = -DV X S Z + ^S x - ^^Vyn, (14) 



Jyx = —DVyS x — OriVn, 



J u 



(15) 



(16) 
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and 



-DVyS z + ^Sy 



-Van. 



(17) 



Now we compare the Eas. (|12[) - (I17|1 with those which 
would be obtained, if we would simply read off the spin- 
current tensor from the di verg ence of the diffusion equa- 
tion, as suggested in Ref.[lOj. Doing so, we note the 
following: while the result for J xy and J yx is in both 
cases the same, the results for the remaining compo- 
nents of the spin-current tensor are different from each 
other (Note, that the fifth term on the lhs of Ea. (|10f> can 
also be written as the divergence of a second rank ten- 
sor). There are two reasons for the differences. Firstly, 
the third term in the Eas. l|14|l and (|17|) has just the 
form of a curl and thus can never be retrieved from 
the diffusion equation. Secondly, the terms linear in 
N in the Eas. ljT^ - ijTTI) differ from those in the diffu- 
sion equation by a factor 1/2. This factor can be traced 
back to the following fact: the second, the third and 
the fourth term in Ea. l)l()|l can be written in the form 
(n • S - DAS - lu s {N x V))j - D llm DD lmk S k , where 

Dilm = OlmVi H eikm e kiz (18) 

h 

is the covariant derivative, the image of the velocity op- 
erator in the diffusion formalism. Thus, the spin diffu- 
sion equation IjlOfl has a similar structure as the time- 
dependent Ginzburg-Landau equation. In Ea. H10[) the 
magnetic field is replaced by the second term on the rhs 
of Ea. (|18l) . Collecting all terms, which are proportional 
to the derivative, means counting the connection form 
twice, a fact which is known to be incorrect from the 
Ginzburg-Landau equation. 

The Eqs.(H2J|-l|n|) differ from those of Ref. in the sec- 
ond term on the rhs of the Eas. i|13fl and IjlSp. This term 
is absent in Ref. 4]. However, it has important conse- 
quences: since the particle density does not vanish at the 
boundary the magnetization can not be zero there. To in- 
vestigate this point further we now focus on the situation 
in the half plane y > 0. Doing so, we ignore the im pac t 
of the spin-orbit interaction on the particle density [l2(. 
First we focus on the situation in equilibrium. In this 
case V x = 0, d/dt = and the particle density n is equal 
to a constant no. Solving the spin-diffusion equation l|l(J|) 
subjected to the boundary condition J y k\ y =o — yields 
Sy — S z — and 

fn , fn . 

S x = tNhq y — exp(-W — y). 



(19) 



Thus, S x is non-vanishing at the boundary. If instead of 
half plane we would have chosen a strip of finite width 
2b wc would have obtained 



S x = -rNn 



IT sinh( 1 



%y) 



instead. 

Now we bias the system. To this end we use the 
method of Ref. 0] , that is instead of switching on a true 
electric field we introduce a concentration gradient. In 
this case the stationary solution to the particle diffusion 
equation is given by n(x) = tiq + n\x. The presence of n\ 
leads to a spin-accumulation, as discussed in Ref. 4] . Due 
to n\ we obtain S x = S z = and S y — Sbuik — —riirN 
in the bulk. Despite this fact, the spin-Hall current van- 
ishes in the absence of a boundary, since the second term 
on the rhs of Ea. (|17fl cancels the third. To investigate the 
situation in a half-plane we impose again the boundary 
condition J y k\y=a = 0. In this case, we obtain 



S z = -SbuikSgnN 
,(1 + A* 2 )(1- 



Re(- 



a 2 )(a* - i) , , m 

A *_ A A W-A^/- y ))(21) 



cosh( 1 



Here A = y— l/2 + i^/7/4, and * is complex conjuga- 
tion. The result for S x agrees with Ea. (|19|l up to a re- 
placement n — > n(x). If instead of the half plane we 
consider a strip of finite width we obtain again a result 
which is antisymmetric with respect to y. 

Here, however, the question arises what the source of 
the spin-Hall effect is, if the spin-Hall current vanishes in 
the bulk. The answer is, that the spin-Hall effect is not 
driven by the spin-Hall current, but by the equilibrium 
magnetization, thus finally by J yx . Once the magnetiza- 
tion exists the field induced rotation (the fourth term in 
Eq.(JT3) turns the magnetization into the z-direction if 
the system is biased. 

The existence of the equilibrium magnetization, how- 
ever, can not be reconciled with time reversal invariance. 
The only way to restore time reversal symmetry is to as- 
sume that the spin-current density has a discontinuity at 
the sample boundary. The discontinuity results from the 
coupling of the particle spin to the particle momentum, 
which is provided by the electron spectrum. Owing to 
the coupling the spins of the particles running toward 
the boundary are mainly aligned in positive ^-direction. 
The particles, which are reflected at the sample bound- 
ary have their spins aligned in the opposite direction. 
Thus, the spins flip when the particles are scattered elas- 
tically at the boundary. This fact has to be taken into 
account in the interpretation of the spin-current tensor 
Q). Whereas the electric current across the boundary is 
a measure for the difference between the number of par- 
ticles running toward the boundary and the number of 
particles reflected at the boundary the spin current ten- 
sor is actually a measure for the sum of both, for the 
number of particles running toward the boundary and the 
number of particles reflected at the boundary. Therefore, 
the spin-current tensor Q does not characterize a flux 
of spins across the boundary in the conventional sense in 
(20) the Rashba model and jumps at the sample boundary. 
The correct boundary condition is S {boundary — in the 
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present situation. In this case there is neither an equi- 
librium magnetization nor a spin-Hall effect in a dirty 
Rashba semiconductor. 

Clearly, at this point the question arises why there 
should be a non-vanishing spin-Hall effect in an electric 
field in the same model in the clean limit, in which the 
energy level splitting due to the Rashba interaction is 
large compared to the disorder energy. The electric field 
itself does not affect the time reversal invariance. The 
boundary condition, however, does in this case. Despite 
this fact energy conservation requires that the spins have 
flipped when the particles have been reflected elastically 
at the sample boundary, at least if the particles are sep- 
arated from the boundary by a distance of the order of 
the mean free path. Consequently, we expect that the 
spin accumulation is confined to a thin layer, which has 
at most a width of the order of the mean free path. 

In summary, we have shown that the assumption 
that the spin-current tensor is continuous at the sample 
boundary leads to results, which arc in contradiction to 
time reversal symmetry. The only way to reconcile our 
results with time reversal symmetry is to assume, that 
the spin-current tensor has a jump at the sample bound- 
ary. The existence of the jump can be traced back to 
the coupling between spin and charge. We would like to 
stress that our results are specific to systems with Rashba 
spin-orbit interaction and do not apply to systems, in 
which the spin-orbit interaction is produced by the scat- 
tering at impurities. The latter differ from the systems 
discussed here in that there is no coupling between spin 



and momentum in the spectrum. 
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A consideration of the impact of the spin-orbit interac- 
tion on the particle density does not effect any of our 
arguments. It leads to a weak accumulation of charge 
carriers at the boundary, which does not affect the final 
result. 



